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We examine the XXZ model with generalized periodic boundary conditions and identify conditions for 
the truncation of the functional fusion relations of the transfer matrix fusion. After the truncation, the 
fusion relations become a closed system of functional equations. The energy spectrum can be obtained 
by solving these equations. We obtain the explicit form of the Hamiltonian eigenvalues for the special 
. case where the anisotropy parameter q 4 = — 1 . 
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^ ; 1 Introduction 

O 1 „ 

The Hamiltonian of the anisotropic Heisenberg-Ising model, or the XXZ model, has the form [jy 

C^: N / q + q- 1 z z \ 

P • HXXZ = ^2 ( a n a n+l + a n a n+l H 7 KK+l J ■ (!) 

43 ■ 

We focus on the case where the module of the anisotropy parameter q is equal to 1, i.e., q = e lr) . We 
^ ' consider the generalized periodic (twisted) boundary conditions 

'x : 

H a N+l 



= q^af (2) 



The spin projection S z = \a\ + ... + \a 7 ^, on the z axis is conserved, i.e., [S z , Hxxz] = 0. As 
established in Q , the energies of chains with different twist parameters are interconnected such that the 
energy spectrum in the chain with the twist parameter (3 from the sector with S z = f3 — 1 contains the 
energy spectrum of the chain with the twist parameter f3 — n from the sector S z = (3 — 1 + n, where n is 
an integer, i.e., 

F (P) _ p(/3-«) 

a S'=0-l ~ rj S*=f3-l+n- 

This properly follows from the quantum group symmetry U q (sl(2)). We say that "one spectrum contains 
the other" because the number of vectors in the sector S z = (3 — 1 + n is smaller than in S z — (3 — 1. 

The XXZ model is connected with the two dimensional classical statistical lattice six vertex model 
(the ice model) ED]. The transfer matrix {^^(u) of the six vertex model commutes with the Hamiltonian, 
$1/2(11), Hxxz] — 0, for all values of the spectral parameter u. This matrix also commutes with itself, 
[ti/2{ u )i ti/2( v )] = 0, and is therefore a generating function of the commuting integrals of motion that are 
in involution. In particular, the Hamiltonian itself can be expressed through the logarithmic derivative 
of t\/2(u) as 

Hxxz = ~ cos(r?) + sm(r))-^- logi 1/2 (u)\ u=0 . (4) 

Formula (^) allows obtain the energy spectrum of the XXZ chain if the transfer matrix eigenvalues are 
known. 
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There exists an infinite family of commuting transfer matrices tj(u), where j — 0, ^, 1, 1, . . .. These 
matrices satisfy the infinite system of recursive functional fusion relations || [|. |[ 

h/2{u-(j + l/2)rj)ij(u) = t (u- (j + 1)77) i 3 _ 1/2 {u + r,/2) 

+ t (u- jrj) i j+1/2 (u-rj/2), (5) 

where to(u) = sin N (u + r//2). In this article, we prove that in the case where the anizotropy parameter 
is the root of unity, q p+1 = —1, the transfer matrix in the representations with the spin j = p/2 has 
zero eigenvalues on the Bethe eigenvectors if the twist parameter j3 and the sector S z — s satisfy the 
conditions 

(-1) 2/3 = {-l) N , \s\ <mm(/3, p+ 1-/3) mod (p+1) (6) 

Once the eigenvalue of t p / 2 {u) vanishes, the infinite system of functional fusion relations (JsJ) is truncated 
and becomes a functional equation with respect to ti/ 2 (u). Solving this functional equation, it is possible 
to obtain the XXZ energy spectrum through formula (Q). We let the symbol V p ^- S denote the set of 
the Bethe eigenvectors on which t p / 2 (u) has zero eigenvalues. In this article, V p ^- S is fully determined in 
the terms of the generators of the quantum algebra U q (sl(2)). We call the XXZ model reduced onto the 
set V p fi- S the reduced XXZ model. The energy spectrum of the reduced XXZ model can be obtained 
by solving the functional equations for the transfer matrix of the siz-vertex model. 

In Sec. 2, we prove the existence of the Bethe vectors on which the functional fusion relations are 
truncated if conditions (g) are satisfied. All such vectors are identified in Sec. 3, where we define V Pl p-s 
in terms of the generators of the quantum algebra U q (sl(2)). In Appendix A, we derive the formula for 
the transfer matrix t p /2{u). The transfer matrix eigenvalues and the energy spectrum of the reduced 
XXZ model for the case where the anisotropy parameter q A is equal to —1 are explicitly calculated in 
Appendix B. 



2 Zero eigenvalues 

We prove the existence of zero eigenvalues of the transfer matrix t p / 2 (u) assuming that q p+1 = -1 and 
conditions (||) are satisfied. The transfer matrix of the siz-vertex model corresponding to the XXZ model 
Hamiltonian with twisted boundary conditions has the form 

t 1/2 (u) = tr (q-P°*R N (u) . . . , (7) 

where the i?-matrix is 

p r„.\ ( sin(« + 1(1 + ct*)t/) singer" \ , . 

Kn[ ] \ MriX + 1(1 - Or,) ) ■ { ' 

The well-known procedure of the algebraic Bethe anzatz (quantum inverse scattering method) Q leads 
to the transfer matrix eigenvalues expressed through the roots {v m } of the Bethe equations, 

n , , tt sin(u — Vm — v) 

tl/a( «) = q -%(u + v/2)l[ i (u _ m Vm ) ^ 



rn—1 



+ q p t {u-r)/2) II — — , 

sm(u — v m ) 

where n — — s is the number of the flipped spins in the sector S z — s, and to(u) is the transfer matrix 
in the representation with the spin j = 0. The Bethe equations are equivalent to the conditions that the 
transfer matrix eigenvalues determided by Eq. (^) have no poles in the complex plane of the variable u. 
In the variables 

n 

Q(u) = JJ sm(u-v m ). (10) 

m— 1 

Eq. (^) is equivalent to the Baxter T-Q-equation 

t 1/2 {u)Q(u) = q- p t (u + rj/2)Q{u - V ) + q rJ t (u - t]/2)Q{u + 77). (11) 

As shown in Q , it is convenient to interpret this equation as a discrete realization of a second-order differ- 
ential equations M. In additions to the functions Q(u), this equation should then have a second lineary 



independent solution P(u) with the same eigenvalue of the transfer matrix t(u). The transfer matrix 
eigenvalues can be expressed through the eigenvalues of the Baxter operator Q(u) and the operator P(u) 



t 3 (u) = qW-Wf{u-{j-l/2)ri)x (12) 
x {Q{u-{j + l/2) V )P{u+{j + l/2) V )- 
- Q(u+(j+ l/2)n)P(u - (j + l/2) V )] , 

where the function f(u) is quasiperiodic, 

f{u + r l )=q-^f{u). (13) 

This function is introduced for convenience (the function P(u) = f(u)P(u) could be introduced instead). 
Because the functions tj (u) satisfy the same functional fusion relation (|5|) , we conclude that these functions 
are some eigenvalues of the corresponding transfer matrices tj (u) . 

We now consider the case q p+1 = — 1, where p is a natural number. For the transfer matrix in the 
representation with the spin j = p/2, we have 

t p /2(u + ~) = {-lff{u)[P{u + n)Q(u) - P(u)Q(u + tt)]. (14) 



It then follows that if the functions Q(u) and P(u) satisfy the relation 

P(u + tt) _ Q(u + tt) 
P{u) ~ Q{u) ' 



(15) 



then the eigenvalue of the transfer matrix t p /2(u) is zero. 

By analogy with the procedure in |Q, we can express the function P(u) through the function Q(u). 
We use the expansion 

to(u) A(u + v /2) _ A(u-r)/2) 

Q(u + V /2)Q(u- v /2) a[U) + q Q(u + v /2) 9 Q(u-r//2)' 1 °> 

where R(u) is a trigonometric polynomial of the order N — 2n, uniquely determined through the known 
trigonometric polynomials to(u) and Q(u), and the order of the polinomial A(u) is smaller than n. 
Further, we use another expansion, 

R(u) = q F(u + ?7/2) - q-PF{u - T]/2). (17) 

It is now easy to verify that the function P(u) is determined by the expression 

P(u) = -]-(Q(u)F(u) + A(u)). (18) 

This expression implies the following constraint imposed on the values of the twist parameter (3 in the 
boundary conditions of the XXZ model: 

("1) 2/3 = (19) 

This means that if the length of the spin chain of atoms is even, then the number must be integer, and 
if the length is odd, the number must be half-integer. No truncation of the functional relations occurs 
for other values of (3. Representation ( |l7| ) of the function R(u) and F(u) in terms of the functions F(u) 
in which F(u) behaves exactly as R(u) under the argument shift u — > u + 7r exists if the inequality 

f} Bm(n^±l) f 0, (20) 
J -- L p + 1 

m— — s 

holds, becase this product appears in the denominator if F(u) is expanded in the Fourier series. The 
variable m takes either integer or half-integer values from — s to +s. Depending on whether the chain 
length is even or odd, the numbers f3 and s are either both integer or both half-integer, and their sum or 
difference is therefore always integer. Condition (|2^) is satisfied if 

\s\ < min(/3, p+ 1 - (3) mod (p+ 1). (21) 

Hence, analyzing the second lineary independent solution of the Baxter equation, we conclude that for 
the anisotropy parameter equal to the root of unity, q p+1 = — 1, the transfer matrix with the spin j = p/2 
in the sector S z = s has zero eigenvalues if conditions (fl9) and (pL) are satisfied. 



3 Algebraic structure 



Above, we proved the existence of vectors such that functional relations (||) restricted to these vectors are 
truncated and have form of a closed system of equations for the eigenvalues of the transfer matrix ti/ 2 (u). 
In this section, we find all vectors from the space of states of the model that satisfy this condition. 

As mentioned above, the connection between the energies of the chains with different twist parame- 
ters (||) exists because of the presence of the quantum group symmetry U q (sl(2)). The generators of this 
symmetry are the operators S z , A and A*, where 



N 



x = ?'K+-+<-i) (T + ? -lK +1 +-+^) (22) 

n— 1 

It is easy to verify that XP+ 1 = for qP +1 = -1. 

The operator A acting on a vector from the sector S z = s transforms this vector to a vector from 
sector S z = s + 1: [S z , A] = A. Acting with this operator on the Hamiltonian, we obtain 

X Hf xz A = (. . .)(1 - q^-^). (23) 

In this formula and in Eq. ( |32| ) below, the dots stand for a finite number of operators whose explicit form 
is cumbersome and irrelevant to this consideration. Formula ( |2"g| ) proves relation (||). A similar relation 
holds for the transfer matrix. 

The six vertex model transfer matrix related to the XXZ Hamiltonian is expressed through the 
matrix elements of the monodromy matrix L(u) = Rn(u) . . . Ri(u), 

2j+l 

t ] {u)=Y J 1- 2P(3+1 - n) Ll(u). (24) 

n=l 

The i?-matrix in the representation i (8 j on lattice site with the number a has the form 

TtUA-f sin ( u + (I +H a )r}) sm(r])F a \ , . 

a[ ' ~ \ Mv)E a sm(u + (| - H a )r,) J • 

The operators E, F, and H have the matrix elements 

Kj(fi)™ = (j + l-n)5 m , n , m, 7i = l, 2, 2j + l. ^ 

ltj{E)™ = UJ m 8 m ,n-l, Kj{F)n = w " S m -l,n, 



where w„ = yl [n] q [2j + 1 — n] q In this formula, we use the standard notation [x] q = (q x — q x )/{q — 
q _1 ). In particular, the operator H for j — 1/2 is ^cr 2 , and the operators E and F are respectively a + 
and <7~. (i.e., the generators of the quantum algebra U q (sl(2)): [H,E] — E, [H,F] — —F, [E,F] = 

m q ). 

Using the Yang-Baxter equation 

(R i9j )Z(u - v){L k ) b c {u){L 3 )™{v) = (L 3 r m {v)(L k )l{u){R k ^(u - v), (27) 

as in ||], we obtain the commutation relations between the matrix elements of the monodromy matrix 
and the operators S z and A, 

q s *Ll = q n ~ k Llq s \ (28) 
XLl{u) = q 2 ^- n - k Ll{u)X 

+ Wfc _ 1 e + »^ +1 -"^ 1 ( U )g- 5= 
- u n e+ i W+ 2 - k L n k +1 (u)q+ s *. 

To obtain commutation relations d28|), we consider the monodromy matrix (Li/2)k( u ) as u — > — iao. 
The matrix element (Li/2)\{ u ) is proportional to the operator A. We then pass to the limit u —* —ioo 
in Yang-Baxter equation (p7f). Commutation relations ( p8| ) allow expressing the transfer matrix itself 
through the operator A. The formula for the transfer matrix with the spin j = p/2 in the sector S z = s 
is (see Appendix A) 

v 

t P / 2 (u) = xr M rk X k , (29) 

r,k=0 



where 



M, 



J-rk 



P 

q(2{fi-s)+r-k)nQn(jp-n 

n=0 

p— i — 

^ ^"i:P-»--fc(ip/2)( u )p-r-fc-m^ P 
m=0 



(30) 



u) = Ylk=i ^ e numbers C% are the g-binomial coefficients 



a 



N 



[N] q [N-l) q ...[N+l 



WW- 



J9 ! 



] g ![JV-n] 



and the variables \ m ,p~r-k are determined by the formula 

m 

\ imu p/2 J rl+s+{m > — m) (m' + s — 3p/2 — 1) rvm I I 



1 = 1 



As can be seen from formula the transfer matrix t p / 2 {u) is expressed through of monomials of the 



'rk^ 



We now investigate the numbers r and k for which the operator M r ^ is nonzero in the case where 



■1. The operator is proportional to the sum of the g-binomial coefficients 



fr 



(2l+r-k) 



(31) 



n=0 



where I = (3 — s. Unfortunately, we failed to obtain an expression for the coefficients in a closed 
form. We therefore analyzed sum ([n]) numerically. We verified the hypothesis that for q p+1 = — 1, the 
coefficients f r k are nonzero only in the closed domain r + k < p and (r > p — I or k > I — 1). 



p-r 











frk = 


r + 








V- /rfc = 
Mill, 


1 1 1 


N l 1 » 



P k 



Hence, the transfer matrix t p i^{u) can be expressed as 

t p/2 {u)=X^-\...) + (...)X l 
and therefore vanishes on the cohomologies 

V p ,i = KeiX l /lmX p+1 - 1 . 



(32) 
(33) 



This means that t p /2(u) has zero eigenvalues on the eigenvectors v form the sector S z = s that are 
annihilated by the operator X, 

X 1 v = 0, (34) 
and cannot be expressed in the form X p+1 ~ l x, where \ is again vector, 



v ? X p+1 -^-^x- 



(35) 



We have thus identified all vectors on vich the transfer matrix t p /2{u) has zero eigenvalues and, as 
mentioned above, recursive system (^|) takes the form of a closed system of equations on the eigenvalues 
of the transfer matrix ti/ 2 {u). 

An example illustrating how the truncated functional relations can be used is given in Appendix B. 
Specifically, the eigenvalues of the transfer matrix and the Hamiltonian arc obtained there for the case q A — 
-1. 



4 Appendix A 



Formula (|29|) can be directly derived by using commutation relations (|28|), but this way is rather weari- 
some. The difficulties are connected with the action of the operator X on L2(u) resulting in the appear- 
ance of two terms in the right-hand side, namely, L^_ 1 (u) and L^ +1 (u). The calculation, however, can 
be simplified. We note that the commutation relations for X and L?(u) are formally equivalent to the 
commutation relations for X and the auxiliary operator %p n {u)q~ nS ip m (u)q~ mS , where the auxiliary 
objects if) commute with X according to the rule 

XipJ - q j+1 - 2n ipJX = -w n e iu q^ +1 (36) 

q S ^i = q ]+1 - n ^ n q SZ 
With this definition of the auxiliary objects ip, the correspondence rule 

(£,-)£(«) ~ q-" 8 ' i>{{u)q- k3 \ (37) 

can be formulated. The symbol ~ means that the commutation relations for the variables X and F(Lu), 
where F is an arbitrary function, are exactly the same as those for X and F(ip n (u)q~ nS ipk{u)q~ ks ). 
The commutation relations for the latter variables, however, are easier to obtain because commutation 
relations ([36]) are simpler than relations (pSj) . Hence, the symbol ~ in the Eq. ([37]) should be interpretted 
as "transforms as" . Hereinafter, we use the equality sign instead of the symbol "~" for simplicity. 
Relations (36) can be rewritten as 

r 3 = -K-ig)-^- 1 '"^™- 1 -^^ 1 - 2 ^- 1 ^"" 1 ^) ( 38 ) 

< = ^ n+ x)- l e- iu {X^ n+1 -q^-^^ n+1 X). 

These simple recursive relations allow expressing ip n through ip 1 and ip n through ip2j+i- The solution of 
thes equations have the forms 

71—1 

W = a n Y / {-l) r q rU+1 - n) C r r l - 1 X n - 1 - r ^ 1 X r (39) 

r=0 

2j+l-n 

^ = b n (-i) fe 9- fe(j+1+n) cf +i -"x^+ i -"- fc v^ + ix fc , 

k=0 

where the coefficients are 

n-i 



Hi-e-^q-^- 1 ) (40) 



k=l 

2j+l-n 



II (e — , 



"2j+l-fc) 
k=l 



The transfer matrix is 



2j+l 

tj(u) = ^Q' 2 ^ 1 ' 71 ^ 1 ^) (41) 
n=l 
2j+l 

\^ q- 2 W +1 - n ^(u)q- nSZ iPi(u)q- nS * . 



n=l 



According to formulas (B9h, we obtain 



t p/2 (u) = w -i e -P>-2 S - P ( 2/3+ i)/ 2 ^^^ ( _ 1)IV p, / 



-l p -'P%-2s-p(2/J+l)/2 ^ ^ ^(-1 r„-P r /2-W2+2)fc 
n=0 r=0 k=0 

X q nW-s)+r-k) C n C P-n x r ^ X p - r ~ k ^ p+l ) X k . (42) 



It remains to remove the auxiliary operators ip 1 and ip p +i that participate in the formula in the form 
ip 1 X p ~ r ~ k ip p+ i. To do this, we commute the operators X p ~ r ~ k and ip p +x. Using (|3^), it is easy to prove 



by the induction that 

M 

^ A X M ll> j B = e imUq(M-m)(j+l+m,-2B)+A(] + l-B+m) 



]>b < L(u)i_ m q( A + B -™ )S*X M -™ . (43) 



\i=i 



Using this formula and taking the relation [S z , X] = X into account, we obtain expression ( p9| ) for t p / 2 (u). 
Only X and L participate in the final formula; the auxiliary objects ip are no longer needed. 



5 Appendix B 

We obtain the eigenvalues of the transfer matrix and the Hamiltonian for q 4 = —1. We have 

*3/2(u) = 0, (44) 
t w (u) = -(-l)*+' 3 -%( U + 7r/2). 

To simplify the calculation, we introduce the function Sn(u) = (— 2) A '/ 2 i 1 / 2 (u — rj/2). The dunctional 
relations become 

S N (u + l)S N (u - J) = cos^u) - sin^pu). (45) 

8 8 

Equation ( [45| ) contains the numbers iV and I = f3 — s. We obtain four different equations, depending on 
whether these numbers are even or odd. Correspondigly, we obtain four different solutions: 



A I 



Sl${u) = n( 1 + (- 1 )' imcos ^^^) cos ( 4u ))' ( 46 ) 

m—l 

_ M 2 

5|m } +2 («) = sin(2 U )^n( 1 + (- 1 )" mcos ^T^q^) cos ( 4u ))' 

771=1 

j-j- j, e (_ 1 )«m_gm_ gin ^ 2w + — ) - ie _(_1) " m 4i^T s in(2u - ^)), 



2M-1 



m—l 



^ + i(«)=e- tu (-i)'~x 
]f[ (e^ _1 )™ m 3^TT sin(2u + ^) - ie - ^ -1 ^""" sin(2tt - ^)), 



2i\/ 



m— 1 

In this formula, the set {n m } is any set of numbers equal to either or 1, for example, {0, 1, 0, . . . , 0, 1}. 
We interpret them as the fermion occupation numbers. The functions written above satisfy the equation 

S iM+ s(u + ^)S iM+ s(u - |) = cos 4M + 5 (2u) + r s S m 4M+s (2u). (47) 

Hence, we obtained the solutions only for the case where f3 — s is an odd number. For the energy 
eigenvalues, we obtain 

s n \ AM 4 , , 2m— 1, , , 



2^2 V2 ^ 

v v m—l 



= + l-2^(-l) ^cos(7r 



2\/2 v 7 ^ \ ^ 2M+1' 



m—l 
2M-1 



w _ 4M-1 1 ^ -2i(-l)"-m 

/„> 4M + 1 if ^ ,-2i(-l)" TO 7rmA 

K?/, , = ^ = \i + 2i > exp( i — '- . 

4M+1 2V2 y/2 \ ^ 1 4M + 1 ^ 

Although the imaginary unit explicitly participates in the last two formulas, the energy E is real in all 
the cases because the imaginary parts cancel each other in the sum. For spin chains whose length are 
even, N = 4M and N = AM + 2, we obtain 2 M energy levels. For spin chains whose length are odd, 
N = AM - 1 and N = AM + 1, we respectively obtain 2 2M_1 and 2 2M eigevalues. 
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